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The fundamental symmetries of .jV = 1 super Yang-Mills (SYM) theory have two opposite 
thermodynamical behaviors. Scale and chiral invariance are expected to be restored at sufficient 
high temperatures, but on the contrary supersymmetry itself is explicitly broken unless the temper¬ 
ature is zero [1]. Supersymmetry is therefore expected to play only a little role in the thermody¬ 
namics of SYM models and results can be qualitatively extended to strongly interacting theories 
without supersymmetry. 

The jV = 1 super Yang-Mills theory shows confinement at low temperatures and chiral sym¬ 
metry is spontaneously broken, like in QCD [2, 3, 4]. However in jY = 1 SYM both phenomena 
are realized exactly and define two distinct phases separated by a critical point: the relation between 
deconfinement and chiral symmetry restoration can be studied exactly, because in this theory there 
are exact order parameters for the transitions, see the following sections. In QCD chiral and decon¬ 
finement phase transition occur approximately at the same temperature, but it is unclear whether 
this fact is a coincidence due to the lack of an exact order parameter [5, 6]. Numerical simulations 
of jY = 1 SYM can provide thus a better understanding of the thermodynamics of strong interac¬ 
tions and can test the correctness of the effective models that have been proposed so far, see for 
instance Ref. [7, 8, 9]. 

In this contribution we present the first lattice results about the phase diagram of .jY = 1 
supersymmetric Yang-Mills theory with gauge group SU(2) at finite temperatures. The detailed 
discussion regarding the main simulations have been published in Ref. [10]. 

1. The model 

The jY = 1 supersymmetric Yang-Mills (SYM) theory is a SUSY model that describes strong 
interactions between gluons and “gluinos”. The gluino is a spin-'/z Majorana fermion in the adjoint 
representation of the gauge group and it is the superpartner of the gluon. The action on-shell S 
reads in the continuum 

S = Jd\(^-^Tr{F^^Fn + ^X{x){iIp)X{x)^ , (1.1) 

and it is equal to the action of one flavor QCD if the gluino A (x) is replaced by a Dirac fermion 
field in the fundamental representation. 

In our simulations the gauge part of S is discretized on the lattice with the Symanzik tree- 
level improved action Sws, while 0 is discretized using the Dirac-Wilson operator Dw- A gluino 
mass term is introduced to fine tune the lattice theory to a consistent continuum limit where the 
renormalized gluino mass vanishes and supersymmetry is dynamically restored [11, 12]. 

The gluino fermion field can be integrated out in the partition function as in QCD but the result 
is the Pfaffian of the Dirac-Wilson operator 

.r = jDU^{x)Pf{CDw)exp{-l5Sws}- (1-2) 

The Pfaffian can be positive or negative unlike the determinant of this theory and therefore the 
integrand above cannot be interpreted as a probability weight. The Pfaffian is equal up to its sign 
to the square root of the determinant: configurations are generated by the Hybrid Monte Carlo 
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Figure 1: Expectation value of the Polyakov loop and its susceptibility as a function of (C = 2 {ml+ 4 ) ^ 

lattice 12^ x 4 with j3 = 1.65. 


algorithm accordingly to the distribution 

J"' = jDU^{x)y/det{CDw)exp{-l5Sws}, (1-3) 

while the sign of Pf(CDw) is inserted as reweighting factor in Monte Carlo simulations. A sign 
problem will appear if (sign(Pf(CDw))) = 0. The sign of the Pfaffian has been measured and it is 
normally always positive if the gluino mass is sufficiently far from the chiral limit. 

2. The deconfinement phase transition 

The Polyakov loop is an order parameter for the deconfinement phase transition. Center sym¬ 
metry is preserved in the fermion action of our model, unlike in QCD, since the Dirac-Wilson 
operator acts on the gluino field in fhe adjoin! represenfafion of fhe gauge group. The Polyakov 
loop is fhus an exacf order parameter for any value of fhe gluino mass. 

The bare paramefer space has been explored af finite femperafures varying fhe bare gluino mass 
niQ af several fixed values of fhe gauge coupling g fo search fhe poinf where fhe deconfinemenl phase 
fransifion occurs. For example, in Fig. 1(a) fhe expecfafion value of fhe modulus of fhe Polyakov 
loop (|Pl|) is shown as a funclion of fc = 1/(2mo + 8) for a fixed j3 = {2Nc )= 1.65 on a laffice 
12^ X 4; (|Pl|) acquires a non-zero value af fC Ri 0.155. The susceptibilify of fhe Polyakov loop 
Xp is measured fo precisely idenfify fhe poinf where fhe center symmefry is sponfaneously broken. 
The peak of Xp is a clear signal of fhe deconfinemenf phase fransifion, see for example Fig. 1(b). 

The order of the phase transition can be estimated from the maximal value of the susceptibility 
Zmax(T) considered as a function of the spatial volume V. The exponent v in the expected scaling 

XmaxjVl) _ fYl\ r? If 

XmM~\vJ ^ 

is equal to 1 for a first order phase transition, to 0 for a crossover and equal to 0.654(7) for a 
second order phase transition in the universality class of the thr'ee dimensional Z 2 spin Ising model 
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(a) Finite size scaling 


(b) 


Figure 2: a) Scaling of the value of the susceptibility of the Polyakov loop at the peak for a lattice N'^ x Nt, 
with = 4 and A(s = {8,12,16}, j3 = 1.65 and K = 0. 16. The colored lines are obtained from formula (2.1), 
using Vi = 8^ for the extrapolation to the other volumes V 2 = 12^ and V 2 = 16^. The purple line marks the 
expected scaling for a first order phase transition, the red line a second order phase transition with x = 0.654 
and the green line a crossover, b) Extrapolation to the chiral limit of the critical deconfinement temperature. 


[13]. As shown in Fig. 2(a), our results are in good agreement with a second order phase transition, 
meaning that the order of the phase transition is unchanged with respect to pure Yang-Mills theory 
by the gluino interactions in the region explored. In particular, the Svetitsky-Yaffe conjecture [14] 
seems to hold also in our model. 

The peak of the susceptibility in Fig. 1(b) identifies the deconfinement phase transition for 
a given = 1 SYM theory softly broken by a specific bare gluino mass mo ^ —0.875. Many 
differenf simulations have been done for differenl laffice sizes and differenf gauge couplings. The 
peak of Xp defines a critical line of couplings (/3c, Kc) in fhe bare paramefer space. This line can be 
used fo exfrapolafe fhe critical femperafure fo fhe supersymmefric limif. The crifical femperafure is 
defined as 

' ( 2 . 2 ) 


Tc = 


N,a{l3c 


and fhe laffice spacing a{P) is sef af zero femperafure in ferms of fhe scale wq in a mass independenf 
renormalization scheme [15]. A proper definifion of fhe gluino mass is required fo exfrapolafe Tc 
fo fhe chiral limif, i.e. fo fhe poinf where fhe renormalized gluino mass mr vanishes. The adjoin! 
pion mass nia-n, defined in a parfially quenched approach, is proporfional fo fhe square roof of niR, 
and if provides a precise definifion of fhe chiral limif [16]. We sef Iherefore niR = (wonia-Ttf', up 
fo an inessential mulfiplicafive facfor. 

The final exfrapolafion is presenfed in Fig. 2(b). The crifical femperafure decreases approach¬ 
ing fhe chiral limif. The linear til 


(worc)(mfi) = 0.0190(22)m« +0.2432(45), (2.3) 

evaluafed af niR = 0, defines fhe crifical femperafure for fhe = 1 super Yang-Mill fheory 

wo7}|„,,=o = 0.2432(45), (2.4) 

up fo lattice discrefizafion errors. The lattice acfion has been sampled also using one-level sfouf 
smeared links in fhe definifion of fhe Dirac-Wilson operafor, in order fo roughly confirm fhe main 
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Figure 3: Expected phase diagram for chiral symmetry breaking of o/k = 1 SYM with two colors. 

conclusions obtained. No differences have been found in the qualitative behavior; the final result 

wo7;U,=o = 0.242(37) (2.5) 

is compatible with the above extrapolation. We have computed the critical temperature for pure 
SU(2) gauge theory in terms of wo, to have a comparison and an estimation of the effects of the 
gluino interactions on Tc. The result 


woTc|m«=oc = 0.2941(13) 


( 2 . 6 ) 


leads to the ratio 

, =0.826(18). 
rc(pure Yang-Mdls) 

dependent in principle on the common observable chosen to set the scale (in this case wq). 


(2.7) 


3. Chiral symmetry restoration 

Chiral symmetry is generated by the Ua(1) group of transformations 

A-^expjiaysjA, (3.1) 

and it is broken down to the discrete subgroup Z 2 n^ due to the axial anomaly. The group Z 2 a?^ 
is further spontaneously broken down to Z 2 by a non-vanishing expectation value of the gluino 
condensate (AA), the full pattern is 

Ua (1)-> Z 2 ;v,-^ Z 2 . (3.2) 

anomaly (XX) fO 

At zero temperature, the gluino condensate jumps with a first order phase transition crossing the 
critical line defined by fhe condifion = 0 [2, 17, 18]. The gluino mass plays Iherefore a similar 
role fo fhe exfernal magnefic field B in fhe Ising model, in parficular fhe gluino condensafe is in 
correspondence wifh fhe sponfaneous magnefizafion. If fhe analogy wifh fhe Ising model is correcf, 
a second order phase fransifion is expected in fhe supersymmefric fheory for fhe gauge group SU(2), 
occuiTing af a crifical femperafure where fhe Z 4 chiral symmefry is resfored, see Fig. 3. 
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Figure 4: Lattice 12^ x Nt, j3 = 1.7 and k — 0.192; a) susceptibility of the chiral condensate and b) Polyakov 
loop susceptibility. Note that the two peaks roughly coincide. The red line marks the deconhnement tem¬ 
perature obtained from the simulations discussed in the previous section. 


There are in principle no known relations between and the deconfinement temperature 
T^ discussed in the section before. An interesting possibility arises when both phase transitions 
coincide, but it could be as well that chiral and center symmetry are broken or restored at two 
different critical temperatures. 

We have done some preliminary simulations to address this question. The bare chiral con¬ 
densate is simply the expectation value of the trace of the inverse of the Dirac-Wilson operator, 
but it requires an additive and multiplicative renormalization. The subtracted chiral condensate is 
computed to get rid of the additive renormalization constant, while we avoid the calculation of the 
multiplicative renormalization factor by fixing fhe gauge coupling fo j8 = 1.7. The femperafure is 
changed by varying Nf. The laffice chosen has a volume 12^ x Nt, wifh Nt in fhe range 4-12. The 
bare gluino mass is fixed fo fc = 0.192; fhe adjoinf pion mass is arua-n = 0.388(9) and fhe value 
of fhe scale is wp/a = 2.070(38). 

We have compufed fhe chiral suscepfibilify Xc to idenfify as before fhe critical poinf by ifs 
peak, see Fig. 4(a). A direcf comparison wifh fhe Polyakov loop suscepfibilify seems fo supporf 
fhe possibilify fhaf chiral symmefry resforafion and deconfinemenf occur af abouf fhe same fem¬ 
perafure, see Fig. 4(b). We have done some furfher simulafions closer fo fhe crifical poinf and we 
have been able fo see no evidence for chiral symmetry breaking up to a temperature higher than 
approximately twice the deconfinement one. 

4. Conclusions 

We have presented the first numerical investigation of the phase diagram of the = 1 super 
Yang-Mills theory at finite temperatures. We have found that the deconfinement phase transition 
has a good signal for all values of the gluino mass and that the finite volume scaling seems to be 
in agreement with a second order phase transition in the universality class of the three dimensional 
Z 2 Ising model. 

Our preliminary investigations support the possibility that chiral restoration occurs roughly at 
the same temperature of the deconfinement phase transition. From this perspective = 1 SYM 
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appears to be similar to QCD. Further studies closer to the continuum and to the massless limit are 
required to fully support this scenario. 

A possible direction of future research is the study of jY = 1 SYM compactified with periodic 
boundary conditions applied to fermions; first investigations have been presented in Ref. [19]. 
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